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Abstract. We analyzed all published frequencies, i/i and 1^2, of the twin kilohertz quasi-periodic oscillations (kHz 
QPOs) in bright neutron star low-mass X-ray binaries. The two frequencies are well correlated but, contrary to 
recent suggestions, the frequency-frequency correlation is significantly different from a V2 — (3/2) i/i relation. To 
check whether, although not following the the 3/2 relation, the QPO frequencies cluster around a region where 
~ 3/2, we re-analyzed the Sco X-1 data that were used to report that ratio and show that, because the 
distribution of ratios of linearly correlated measurements is intrinsically biased, although the significance of the 
clustering around 1^2/1^1 ~ 3/2 previously reported in the case of Sco X-1 is formally correct, it does not provide 
any useful information about a possible underlying resonance mechanism in this source. Using the same data, 
we then show that the (unbiased) distribution of QPO frequencies is consistent with a uniform distribution at 
a 2.4a level. To investigate this further, we analyzed a larger data set of Sco X-1 and four other sources, 4U 
1608-52, 4U 1636-53, 4U 1728-34 and 4U 1820-30. We find that for all five sources the distribution of the kHz 
QPO frequencies is not uniform and has multiple peaks, which have no analogy in the distribution of points in the 
spectral color-color diagrams of these sources. Finally, we demonstrate that a simple random walk of the QPO 
frequencies can reproduce qualitatively the observed distributions in frequency and frequency ratio. This result 
weakens the support for resonance models of kHz QPOs in neutron stars. 
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1. Introduction 

, The launch of the Rossi X-Ray Timing Explorer (RXTE) 
' satellite has led to the discovery of high-frequency quasi- 
. periodic oscillations (QPOs) in bright low mass X-ray bi- 
' naries (LMXBs) containing neutron stars (see van der Klis 
2000 for a review). These QPOs provide a probe into the 
accretion flow around neutron stars very close to the com- 
pact object, where effects of general relativity might be 
observable. They often appear in pairs at frequencies of a 
few hundred Hz to more than 1 kHz, from which the name 
kHz QPOs was derived. Their frequencies follow rather 
tight correlations with other timing features of the X-ray 
emission (see Ford & van der Klis 1998; Psaltis et al. 1999; 
Belloni et al. 2002). 

There is currently no consensus as to the origin of these 
QPOs, nor on what physical parameters determine their 
frequencies, which have been identified with various char- 
acteristic frequencies in the inner accretion flow (see e.g. 
Stella & Vietri 1999; Osherovich & Titarchuk 1999; Lamb 
& Miller 2003). 
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Recently, pairs of 30-450 Hz QPOs have been ob- 
served from a few black-hole candidates and their fre- 
quencies tend to appear in certain ratios (3:2, 5:3, see e.g. 
Strohmayer 2001a, b; Miller et al. 2001), although there are 
exceptions (see e.g. Homan et al. 2001). Unfortunately, the 
detection of such QPOs in black-hole candidates is rather 
rare, which makes it difhcult to assess the significance of 
this clustering around definite ratios. Abramowicz et al. 
(2003) reported that the ratio of the frequencies in the 
kHz QPOs from the brightest LMXB in the sky, Sco X-1, 
tend to cluster around a value of 1.5, which they interpret 
as evidence for a 3:2 resonance. A mathematical approach 
to the resonance model was presented by Rebusco (2004) 
in order to explain the Sco X-1 results; the discrepancies 
of the data with a pure 3:2 ratio were attributed to the 
action of an additional ad-hoc force. 

In this paper, we present the results of an analysis of all 
published values of kHz QPO frequencies from neutron- 
star systems, including those from Sco X-1. We also ex- 
tracted and analyzed long series of QPO frequencies from 
a few selected systems: Sco X-1, 4U 1608-52, 4U 1636-53, 
4U 1728-34 and 4U 1820-30. From these data, we con- 
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Table 1. List of sources considered in this work, with range of 
observed vi frequency and references 



Source 


Frequency 


References 


name 


range (Hz) 




Z sources 


GX 17+2 


475-830 


[1] 


Sco X-1 


565-829 


[2] 


GX 5-1 


156-627 


[3] 


GX 340+0 


197-565 


[4] 


Cyg X-2 


532 


[5] 


GX 349-2 


712-715 


[6,7] 


Atoll sources 


4U 1728-34 


308-876 


[8,9] 


4U 0614+09 


153-823 


[10] 


4U 1705-44 


776 


[11] 


KS 1731-260 


903 


[12] 


4U 1735-44 


641-726 


[13] 


4U 1608-52 


550 


[14] 


4U 1636-53 


865-920 


[15,16] 


4U 1820-30 


790 


[17] 


4U 1915-05 


224-707 


[18] 


XTE J2123-058 849-871 


[19] 



[1] Homan et al. (2002); [2] van der Klis et al. (1997); [3] 

Jonker et al. (2002a); [4] Jonker et al. (2000); [5] Wijnands 
et al. (1998); [6] Zhang et al. (1998); [7] O'Neill et al. (2002); 
[8] Migliari et al. (2003); [9] Di Salvo ot al. (2001); [10] van 
Straaten et al. (2000); [11] Ford et al. (1998a); [12] Wijnands 
& van der Klis (1997); [13] Ford et al. (1998b); [14] Mendez 
et al. (1998b); [15] Wijnands et al. (1997); [16] Di Salvo ct al. 
(2003); [17] Smale et al. (1997); [18] Boirin et al. (2000); [19] 
Homan et al. (1999). 



elude that the presence of a single fixed ratio between 
the frequencies of these QPOs can be excluded. We do, 
however, find evidence of a clustering of QPO frequencies 
around specific values that might be associated with reso- 
nances. Additional independent observations are needed 
to exclude a random origin of this clustering. 



2. The sample of published kHz QPO frequencies 

We searched the literature for all published instances of 
kHz QPO frequencies. In some cases, no tables are pro- 
vided and the numbers could only be obtained from fig- 
ures. We did extract numbers from such figures only when 
no other frequencies were available from that source. We 
restricted our sample to the detection of two simultane- 
ous kHz QPO peaks and discarded all data obtained via 
a shift-and-add technique (see Mendez et al. 1998a), as 
this procedure does not allow to reconstruct the distribu- 
tion of pairs of frequencies (see Abramowicz et al. 2003). 
The adopted sources with the range of lower kHz QPO 
frequency and references are reported in Table 1. For Sco 
X-1, we used the same data as in Abramowicz at al. (2003). 




Lower kHz QPO frequency (Hz) 



Fig. 1. Plot of upper kHz QPO frequency (v'z) vs. lower kHz 
QPO frequency {vi) for the Atoll and Z sources for which dou- 
ble kHz QPOs have been published (see Table 1). Empty circles 
correspond to Z sources, stars to Atoll sources and filled circles 
to the Sco X-1 data (courtesy of M. van der Klis). The dotted 
line is the best linear fit to the Atoll points, the dot-dashed 
line is the best lineax fit to the Z points (excluding Sco X-1), 
and the thick line represents a fixed 3:2 ratio. Errors of the 
frequencies are typically around a few Hz. 

2.1. Correlation between frequencies 

Figure 1 shows a linear plot of the frequency of the upper 
kHz QPO V2 versus the frequency of the lower kHz QPO 
vi. Stars indicate Atoll sources, circles are for Z sources 
(filled circles for Sco X-1). It is evident that a strong linear 
correlation exists between the two frequencies, although 
slightly different for the two classes of systems. A linear 
fit, y = ax + h, yields a=0.94±0.02, 6=350±15 for Atoll 
sources and a=0.85±0.01, 6=383±8 for Z sources without 
Sco X-1 (see dashed and dotted lines in Figure 1). Errors 
are at a la significance level. These relations, as well as 
the a=0.73±0.01, 6=469±7 found for Sco X-1 alone, are 
clearly not compatible with the two frequencies having a 
constant 3:2 ratio ( a=1.5, 6=0: continuous line in Fig. 1). 

Following Abramowicz et al. (2003), we plot the his- 
togram of the V2/V1 ratio distribution. A clear peak 
around 1.5, as reported by Abramowicz et al. (2003), can 
be seen (see Fig. 2). However, this peak should be inter- 
preted with some caution. The distribution of the ratio of 
two quantities that follow a linear correlation with 6 ^ is 
largely determined by the distribution of the observed val- 
ues of one of the two variables. Consider that for a linear 
relation 

y/x = a + b/x. (1) 

For high values of x (i.e. x > 6), which are the majority 
in Fig. 1, the ratio will be asymptotically close to a, while 
for low values of x the ratio will be higher. To illustrate 
this, we simulated a sample of frequency pairs around the 
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correlation for atoll sources, but with the lower kHz fre- 
quencies uniformly distributed in the 150-900 range. The 
corresponding histogram can be seen in the inset of Fig. 2: 
a rapid increase toward low values of the ratio is clearly 
visible. This increase is the combined effect of the fact 
that the two frequencies are correlated and of the range 
of observed frequencies. 




1.5 2 2.5 3 

V / V 
2 1 

Fig. 2. Distribution of the V'^jv^ values for points in Fig. 1. 
Black corresponds to Atoll sources, gray to Z sources (including 
See X-1). The inset shows the histogram corresponding to a 
simulated uniform distribution of points around the correlation 
for atoll sources (see text). 



2.2. A preferred 3:2 ratio in Sco X-1? 

While the effect described in the previous paragraph can 
(practically) account for the observed peak in Fig. 2. it is 
still possible that the points in Fig. 1 intrinsically cluster 
around the intersection with the 3:2 line, even though their 
own correlation significantly differs from that. This cannot 
be tested for most of the sources in Fig. 1, as Abramowicz 
et al. (2003) point out, as the real distribution in occur- 
rence of the QPO pairs in these sources is hidden by the 
typical selection procedures applied by the respective au- 
thors. It can however be tested from the published Sco 
X-1 data, which have been selected in an unbiased way 
(see Abramowicz et al. 2003). 

First, we produced a histogram of the distribution of 
one of the frequencies, in this case 1^2- This histogram can 
be seen in Fig. 3 and shows a broad excess around 900- 
950 Hz. A fit of the distribution in Fig. 3 with a Gaussian 
model yields a centroid frequency of i'2=930±8 Hz {la 
error). Then, we fitted the vi vs. 1^2 relation (the inverse 
relation of that shown in Fig. 1) with a linear model, ob- 
taining a=1.363±0.019 and 6=-635±18. Using this rela- 
tion, we can assign a corresponding i^i frequency to i/f and 
compute their ratio. The resulting ratio is i?=1.47±0.04. 
Notice however that a Kolmogorov-Smirnov test shows 



that the hypothesis that the 1^2 distribution is constant 
can only be rejected at a 2Aa confidence level. 

81 1 1 1 1 1 

7- 



6- 




(Hz) 



Fig. 3. Distribution of the frequency of the upper kHz QPO 
i/'2 for the Sco X-1 data analyzed by Abramowicz et al. (2003). 

There is therefore some marginal evidence for a clus- 
tering around a frequency corresponding to a ratio of 1.5, 
marginal because of the small number of points in this 
dataset. We note, however, that although we compared 
the distribution of f2 against a uniform distribution, there 
is no reason to expect a constant distribution in QPO fre- 
quency. Also notice that in Sco X-1 twin kHz QPOs are 
observed only when the source is in the upper Normal 
Branch in the color-color diagram (van dcr Klis et al. 
1996) and therefore, by selecting only data with two peaks, 
we introduce a bias toward low frequencies (since QPO 
frequencies correlate with the position of the source in a 
color-color diagram). Indeed, van der Klis et al. (1996) re- 
port single-peak detections above 1090 Hz. In a number of 
sources, a second low-frequency kHz QPO is found when 
the Power Density Spectra are shifted according to the 
high-frequency QPO and then summed, indicating that 
the second QPO might always be present, although at a 
lower rms. 

In order to investigate the actual distribution of kHz 

QPO frequencies in bright LMXB we examined the high- 
frequency QPO detections of four additional sources, as 
well as from an expanded set of data from Sco X-1. 

3. Data analysis 

The selected data and the procedures used are described 
below. 

— Sco X-1: we used the data analyzed by Mendez & van 
der Klis (2000). Power density spectra (PDS) were ac- 
cumulated from data stretches 16 seconds long, then 8 
consecutive PDS were averaged and searched for high- 
frequency (> 250 Hz) peaks. If two peaks were found. 
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Fig. 4. The series of measured QPO centroid frequencies for 
the five sources analyzed. The X axis is simply the PDS num- 
ber, i.e. the frequencies are plotted in sequence. All discrete 
jumps in frequency correspond to gaps in the data. Gaps longer 
than one hour are marked with dotted lines. Typical errors on 
the frequencies are 1 — 5 Hz. 



The resulting frequencies can be seen in Fig. 4, with 
vertical markers indicating gaps longer than one hour. 

For each source, wc repeated the analysis described 
in Section 2.2. The distributions of 1^2 (for Sco X-1) and 
i/i (for the other sources) can be seen in Fig. 5. As all of 
the QPO frequency distributions sccni to be multi-peaked, 
we fitted them with a sum of Gaussians. Two or three 
Gaussians were needed to fit the observed peaks. The re- 
sulting ccntroids Va, i^fj and v-y arc reported in Table 2. 
The multi-Gaussian models used are also plotted in Fig. 
5. 

For each source wc then used the best-fit V2 — 1^1 re- 
lation to compute the frequency ratios Ra, Rf3 and 
corresponding to the Gaussian peaks, with corresponding 
errors. The results are shown in Tab. 2 and will be dis- 
cussed in Sect. 4. Notice that the errors on R(^ , Rp and R-y 
are dominated by the propagated errors on the coefficients 
a and b from the linear fits to the i>2 — vi relations. The er- 
rors on the ratios are therefore correlated, i.e. even though 
the signs of the errors on the ratios are not known, they 
arc most likely the same for all ratios of each individual 
source. 



only the highest one was considered. In all cases, the 
upper kHz QPO was detected (see Mendez & van der 
Klis for more details). The lower kHz peak was then re- 
covered via a shift-and-add technique using narrow (5 
Hz) frequency bins: this procedure allows one to obtain 
an averaged V2 vs. vi relation (see Fig. 1 in Mendez & 
van der KUs 2000, also reported by Miller 2003). The 
output of these procedures are therefore a time series 
of V2 frequencies and a general correlation V2 vs. vi. 
We fitted this correlation with a linear model, obtain- 
ing a=0.786±0.002, &=433±2. 

— 4U 1608-52 and 4U 1728 34: we used the data an- 
alyzed by Mendez et al. (2001). The extraction pro- 
cedure was very similar to that for Sco X-1, with 
the following differences: PDS were accumulated from 
64s data stretches, a variable number of consecu- 
tive PDS were averaged (in all cases less than 20) 
in order to obtain a significant detection, and the 
lower kHz QPO peak was selected, obtaining the V2 
vs. vi relation through shift and add. A linear fit 
gives a=0.810±0.017, 6=418±10 for 4U 1608 52 and 
a=0.897±0.024, 6=420±17 for 4U 1728 34. 

— 4U 1636-53: we used the data analyzed by Di Salvo et 
al. (2003). The extraction procedure was the same as 
for 4U 1608-52 and 4U 1728-34. To derive [a,6], we also 
used the points by Jonker et al. (2002b). The linear fit 
for this source gives a=0.673±0.017, 6=540±14. 

— 4U 1820-30: The extraction procedure was the same 
as for 4U 1608-52 and 4U 1728 34. Preliminary re- 
sults were shown in Mendez (2002). We used all 
RXTE/PCA data from 1996 October 15th to 1999 
March 1st. For this source, we obtained a fit with 
a=0.741±0.039, 6=466±29. 



SCO X-1 4U 1820-30 4U 1608-62 

' 1 100 I — ' ' ' ' 1 I ' ' — 




800 1000 1200 500 600 700 800 900 500 600 700 800 900 



V2(Hz) V, (Hz) V, (Hz) 

4U 1636-53 4U 1728-34 




750 800 850 900 950 600 700 800 900 



V, (Hz) V, (Hz) 

Fig. 5. Distribution of the frequency of the upper kHz QPO 1/2 
for Sco X-1, and of the lower kHz QPO vi for the other sources 
in our sample (see text). For Sco X-1, the black histogram 

represents the frequency distribution of the upper kHz QPO 
in the data by Abramowicz et al (2003). 

As the frequency of QPOs in these systems is known 
to correlate with the position on the color-color diagram 
(see e.g. Di Salvo et al. 2003), it is possible that the non- 
uniformity of the distributions shown in Fig. 5 is the result 
of an analogous non-uniform distribution of points along 
the color-color diagram of the source. In order to check 
for such an effect for each of the sources, except Sco X-1, 
we produced a color-color diagram with each point corre- 
sponding to one of our QPO detections, using the proce- 
dures described in Di Salvo et al. (2003). In the case of Sco 
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Table 2. Parameters of the hnear fits and centroids of the Gaussian distributions fitted to the histograms in Fig. 5 and 
corresponding frequency ratios (see text). All errors are Icr. 



Source 


a 


6 








Ra 


R/s 




name 






(Hz) 


(Hz) 


(Hz) 








Sco X-1 


0.786±0.002 


432.5± 1.5 


936 ± 2 


1025± 3 


1077±3 


1.46 ± 0.01 


1.36 ± 0.01 


1.32±0.01 


4U 1608-52 


0.810±0.017 


418.3±10.3 


594 ± 4 


668± 2 


787±2 


1.51 ± 0.02 


1.44 ± 0.02 


1.34±0.02 


4U 1636-53 


0.673±0.017 


539.8±14.1 


809 ±10 


856± 2 


897±3 


1.34 ± 0.03 


1.30 ± 0.02 


1.27±0.02 


4U 1728-34 


0.897±0.024 


419.7±16.8 


756 ± 2 


874± 3 




1.45 ± 0.03 


1.38 ± 0.03 




4U 1820-30 


0.741±0.039 


466.2±29.0 


601 ± 1 


731± 2 




1.52 ± 0.06 


1.38 ± 0.06 





X-1, this proved not to be feasible as the strong dead-time 
effects caused by the high count rate and the different off- 
set pointings to this source prevented the production of 
a single homogeneous color-color track. We parametrized 
the position on the color-color track following Di Salvo et 
al. (2003) and obtained for each point an Sa value, i.e. 
an index that characterizes the position along the color- 
color track. The distribution of the points in the v — Sa 
plane is shown in Fig. 6, together with the marginal dis- 
tribution along the two axes. It is clear from the figure 
that the multi-peaked distributions of QPO frequency do 
not have a corresponding multi-peaked distribution in Sa- 
For instance, in the case of 4U 1820-30 there is a gap in 
the distribution of QPO frequencies between i^i ~ 600 and 
ui ~ 650 Hz, but there is no gap in the Sa distribution at 
the corresponding values between Sa ~ 0.9 and Sa ~ 1.0 
(see upper left panel in Figure 6). 

Having ascertained that the distribution of QPO fre- 
quencies is not constant, but shows multiple peaks, we 
need to establish whether those peaks are simply the re- 
sult of the random walk of the QPO frequencies. It is 
known that the frequencies of kHz QPOs do not jump 
in frequency, but show a sort of random walk in time (sec 
Fig. 4). In fact, this behavior is not a simple random walk, 
as frequencies are only observed within a certain range. In 
order to test whether such a random movement of frequen- 
cies in time could produce peaks in the frequency distribu- 
tion, we simulated a simple random walk in frequency. We 
produced a series of 10000 frequencies starting at 700 Hz 
and adding a random value between —6 an +6 Hz to the 
previous value in the sequence. The resulting histograms 
for two realization of the random walk are shown in Fig. 
7. Adopting a linear relation between ui and 1^2 as the 
one shown in Table 2 for Sco X-1, for each simulation we 
also produced the distribution of v^jvx. This very sim- 
plified model, which does not include a "re-call" force to 
keep the frequency between fixed bounds, can clearly pro- 
duce multiple peaks in the frequency distribution, which 
for the reasons outlined above results in a sharp peak in 
the distribution of ratios. 

4. Discussion 

We investigated the possibility of a 3:2 resonance in the 
frequencies of the kHz QPOs in bright accreting LMXBs 
from all available data in the literature. We found that the 



strong linear correlation observed between the lower and 
upper kHz QPO frequencies both in Atoll and Z sources 
is not compatible with a single constant ratio, which can 
be excluded with high significance (see Fig. 1). This con- 
clusion is also evident from Rebusco (2004), who suggests 
an ad hoc modification to the resonance model in order to 
explain the observed correlation. We showed that even a 
constant distribution of the frequency of one of the QPOs, 
together with an observed correlation between the two 
QPO frequencies, would result in a peaked histogram of 
frequency ratios. Since resonance models of the QPOs 
predict that the frequency ratios as well as the frequen- 
cies themselves should cluster around specific values, and 
since the distribution of QPO frequencies does not suffer 
from the bias of the distribution of frequency ratios, we 
used the same Sco X-1 data analyzed by Abramowicz et al. 
(2003) to check whether there is evidence for a clustering 
around a particular QPO frequency. We found marginal 
evidence for an excess in the distribution at a frequency 
consistent with a 3:2 ratio, whose significance is difficult to 
assess given the scarcity of available data, and whose na- 
ture might be the result of observational biases. We con- 
clude that although the significance of the narrow peak in 
the distribution of QPO frequency ratios around a value 
of 3/2 reported by Abramowicz et al. (2003) in Sco X- 
1 is formally correct, since the distribution of frequency 
ratios is biased, their result has no power in testing the 
predictions of their resonance models; the peak in the dis- 
tribution of ratios that they find is probably due to the 
distribution of QPO frequencies in Sco X-1 being mostly 
above the intersection with the line of a constant 3:2 ratio. 

The analysis of a sample of high-frequency oscillations 
in five systems (including a more extensive coverage of 
Sco X-1) shows that the V2 — v\ relation is rather simi- 
lar between different systems, and that the distribution of 
the observed frequencies is not flat in the observed range, 
but shows two or three peaks, depending on the source. 
Abramowicz et al. (2003) reported the presence of a sec- 
ond peak, although at low significance, corresponding to 
the six points above ^2=1050 Hz in Figs. 1 and 3. From 
Fig. 5, it is evident that those points correspond to our 
high frequency component v^. Using the measured v\ — 
relations, it is possible to associate a frequency ratio to the 
centroids of these peaks. The lowest peak corresponds to 
a ratio close to 1.5 in four systems, while for 4U 1636-34 
it does not. 
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Fig. 6. Frequency of the lower kHz QPO vi plotted ver- 
sus the position of the source along the color-color diagram, 
parametrized by Sa (see text), for the four atoll sources in our 
sample. The histograms represent the marginal distributions 
of the single parameters. Notice that the v histograms are not 
identical to those in Fig. 5, because not for all points we could 
extract a color. 

There is no a priori reason for such a distribution to 
have a flat shape, since it is known that the fractional rms 
of these oscillations decreases at high and at low frequen- 
cies (see Mendez et al. 2001 and Di Salvo et al. 2003). 
However, the frequency dependence of fractional rms does 
not show any evidence of being multi-peaked, and its max- 
imum does not coincide with any of the maxima of the 
distributions in Fig. 5. 
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Fig. 7. Left panels: Distribution of the frequencies of two ran- 
dom walks in frequency (see text). Right panels: corresponding 
distribution of frequency ratios assuming the vi — V2 relation 
for Sco X-1 (see Table 2). 

We have shown that the multi-peaked nature of the 
distribution of QPO frequencies is not due to a combina- 
tion of the source spending more time on certain positions 



in the color-color diagram, and of a correlation between 
position on the color-color diagram and QPO frequency 
(Fig. 6) . The lack of correspondence of the multi-peaked 
distribution of QPO frequency with the distribution of 
the Sa values contradicts scenarios in which the QPO fre- 
quency and the source spectrum are both determined only 
by the radius of the inner edge of the disk (van der Klis 
2001). Notice that in particular 4U 1820-30, while show- 
ing a good correlation between V\ and Sa , clearly features 
a gap in the distribution of v\ not associated to a corre- 
sponding gap in Sa- 

From the few instances of double high-frequency QPOs 
in black-hole systems there is some evidence for a pre- 
ferred fixed ratio. The data are sparse and rather stable 
frequencies have been observed in these pairs. In the case 
of BH binaries the preferred ratio can also be the result 
of sampling, i.e. the few detections available to date could 
come from the upper region of a correlation similar to 
that in Fig. 1, in which case the observed ratio would not 
have a simple physical interpretation in terms of a reso- 
nance mechanism. For neutron-star systems, where a large 
number of detections of pairs with variable frequencies is 
available, a single fixed ratio can be excluded. The pres- 
ence of an additional external force as in the model by 
Rebusco (2004) (see also Abramowicz et al. 2003b) could 
explain the observed deviation from a fixed 3:2 ratio, at 
the expense of a considerable complication of the model, 
with the addition of an unknown parameter which biases 
the observed values away from a direct measurement of 
a 3:2 ratio. However, the presence of multiple peaks as 
those shown in Fig. 5 would complicate the interpretation 
further, as the presence of other preferred frequencies (or 
ratios) would also have to be explained. 

If the kHz QPOs are due to a resonance mechanism in 
the accretion disk, the multi-peaked distribution of QPO 
frequencies (Fig. 5) suggests that not just the 3/2, but 
more than one resonance is at work. Using the V2 ~ v\ 
relations in Sect. 3, we can assign a frequency ratio to 
each of the peaks of the QPO distributions in Fig. 5. The 
inferred ratios are close to 3/2, 7/5, and 4/3 in the cases 
of Sco X-1 and 4U 1608-52, 4/3, 9/7 and 5/4 in the case 
of 4U 1636-53, and 3/2 and 7/5 in 4U 1728-34 and 4U 
1820-30. However, the errors are sufficiently large that 
more than one identification is possible. 

From the ratios shown in Table 2, it is apparent that 
4U 1636-53 is different from the other sources in that the 
frequency ratios in this source are systematically lower 
than in the other four sources. While for all other sources 
the largest ratio is close to 1.5, the largest ratio in 4U 
1636-53 is 1.34±0.03, significantly smaller than 1.5. From 
the sample of sources in Table 2, 4U 1636-53 is the only 
one that during the RXTE observations used for our anal- 
ysis appeared always in a region of the color-color diagram 
characteristic of sources at high inferred mass accretion 
rates, and the QPOs appeared always at high frequencies 
(Di Salvo et al. 2003). Assuming there is a relation be- 
tween QPO frequency and radius of the inner edge of the 
disk (e.g., Miller et al. 1998), this suggests that during 
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these observations the inner edge of the disk in 4U 1636- 
53 was always relatively close to the neutron star. On the 
other hand, based on their intensity, color-color diagrams, 
and range of QPO frequencies, the other sources appeared 
to cover a larger range of mass accretion rate, which sug- 
gests that in those sources the inner radius of the disk 
covered a larger range than in 4U 1636-53. 

It has already been shown (Stella & Vietri 1999) that 
the identification of the QPOs with the azimuthal and 
the periastron precession frequencies, V0 and v-p = — u-r, 
respectively, reproduces qualitatively the QPO frequency- 
frequency correlations. [We note in passing that the al- 
ternative identification of the QPO frequencies with the 
vertical and radial epicyclic frequencies (e.g., Kluzniak et 
al. 2004) implies frequency-frequency correlations in con- 
tradiction with the observed data. For instance, unless 
another mechanism is at work, this interpretation pre- 
dicts that the frequency difference between the kHz QPOs, 
1^2 — 1^1, increases as the frequencies increase, opposite to 
what is observed in several sources; see e.g., Mendez et 
al. (1998a)]. The ve^Vp identification can also account for 
the range of QPO frequency ratios seen in 4U 1636-53, 
compared to those seen in the other sources. The ratio 
R = vg/u-p increases monotonically with r, and hence de- 
creases as the frequencies increase (notice that, on the 
contrary, the ratio ve/i^r increases toward higher frequen- 
cies). Since in our observations of 4U 1636-53 the kHz 
QPOs never reached values as low as those seen in the 
other sources, the frequency ratio did not reach high val- 
ues either. We expect that when the mass accretion rate 
in 4U 1636-53 decreases, and the kHz QPO move to lower 
frequencies, the ratio of QPO frequencies will increase and 
cluster around 1.5 as in the other sources. 

There is in principle no a priori reason why higher 
frequency ratios such as 2/1, 3/1, 4/1, etc. should not ap- 
pear. With the identification of the observed frequencies 
with those of Stella & Vietri (1999), for a neutron star 
mass in the range 1.4 — 2.2Mq, a frequency ratio R > 2/1 
would only occur when the frequency of the lower kHz 
QPO is <300 - 500 Hz, while a ratio i? > 3/1 requires 
the frequency of the lower kHz QPO to be <150 - 200 
Hz. In none of the observations presented in this paper 
did the lower kHz QPO frequency reach such a low value. 
However, there are some reports in the literature of lower 
kHz QPO being detected in this frequency range. For in- 
stance, the lowest frequencies reached by the kHz QPO 
in the atoll source 4U 0614+09 are vi = 153.4 ± 5.6 
and 1/2 = 449.4 ± 19.5 (van Straaten et al. 2000), such 
that i? = 2.9 ± 0.2. In the Z source GX 5-1, the lowest 
frequencies so far observed for a pair of kHz QPOs are 
ui = 156 ± 23 Hz and U2 = 478 ± 15 Hz (Jonker et al. 
2000), for which i? = 3.1 ± 0.4. 

Finally, we showed that a simple frequency random 
walk could in principle reproduce the observed multi- 
peaked frequency distributions as well as the distribution 
of frequency ratios. Although we have used a very sim- 
ple approach to simulate the drift of QPO frequency with 
time (e.g. the real frequencies do not follow a pure random 



walk, the observations sample many separate intervals of 
the real frequency evolution, etc.), it is clear that a multi- 
peaked distribution of QPO frequencies (or frequency ra- 
tios) could also result from purely random processes. A 
possible test of the idea that the observed QPO-frequency 
distribution in these neutron star systems is indeed pro- 
duced by an underlying resonance mechanism would be to 
accumulate new (independent) samples of frequencies of 
similar size to the ones presented here, and check whether 
the distribution of QPO frequencies in these new samples 
show peaks with centroid frequencies that are consistent 
with those observed here (although some peaks could be 
missing while some new peaks may appear, depending on 
the range of frequencies sampled by the sources during the 
new observations). If this is not observed, we would have 
to conclude that the peaks are due to the random motion 
of the frequency. Since at present we cannot discard this 
possibility, our results weaken the support for resonance 
models of the kHz QPOs in neutron stars. 

5. Summary and conclusions 

We showed that: 

— The relation between the two simultaneous kHz QPO 
frequencies in a large sample of sources, including Sco 
X-1, is significantly different from a, 1/2 = l-Bvi rela- 
tion. 

— Since in all these sources the two QPO frequencies are 
linearly correlated, the distribution of the ratios of the 
QPO frequencies would be peaked even if individually 
the frequencies are distributed more or less uniformly. 
Since resonance models predict that both the ratio of 
frequencies as well as the frequencies themselves shoiild 
cluster around specific values, the distribution of fre- 
quency, and not of frequency ratios, should be used to 
test the predictions of those models. 

— The hypothesis that in the same Sco X-1 data used 
by Abramowicz et al. (2003) the distribution of QPO 
frequencies is significantly different from constant can 
only be rejected at a 2.4(t confidence level. 

— The distribution of QPO frequencies in Sco X-1, 4U 
1608-52, 4U 1636 53, 4U 1728 34, and 4U 1820 30 
is multi-peaked, with the peaks occurring at frequen- 
cies compatible with 1/2/1^1 ratios of 3/2, 4/3, 5/4, 7/5, 
and 9/7, not all ratios appearing in all sources. This is 
considerably different from the claim by Abramowicz 
et al. (2003) that the frequency ratio for Sco X-1 shows 
a narrow peak at 1.5, which we showed in Sect. 2 not 
to be statistically founded. 

— The multi-peaked nature of the frequency distribution 
appears difficult to explain by parametric resonance 
models, since the perturbation mechanism that drives 
the ratio away from the main, 3/2, resonance would 
have to be able, at the same time, to produce the other 
resonances as well. 

— A random walk of the QPO frequencies qualitatively 
reproduces the observed QPO frequency distribution 
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as well as the frequency ratio distribution. This means 
that, until the presence of the same peaks at the same 
locations is confirmed by (large) independent datascts, 
it is premature to speculate further about the presence 
of preferred ratios. 
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